Ferromagnetic ground state of an orbital degenerate electronic model for 
transition-metal oxides: exact solution and physical mechanism 
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We present an exact ground state solution of a one-dimensional electronic model for transition-metal 
oxides in the strong coupling limit. The model contains doubly degenerated orbit for itinerant elec- 
trons and the Hund coupling between the itinerant electrons and localized spins. The ground state is 
proven to be a full ferromagnet for any density of electrons. Our model provides a rigorous example 
for metallic ferromagnetism in narrow band systems. The physical mechanism for ferromagnetism 
and its relevance to high-dimensional systems, like Ri-^X^MnOs, are discussed. Due to the orbital 
degeneracy of itinerant electrons, the superexchange coupling can be ferromagnetic rather than 
antiferromagnetic in the one-band case. 



PACS numbers: 71.27.+a; 75.10-b; 75.30.Et 

The origin of ferromagnetism in narrow band systems 
is a long standing problem and has recently attracted 
renewed attention. Essentially speaking, there are three 
basic routes to ferromagnetism: 1). it is formulated in the 
Hubbard model for an intermediate short-range Coulomb 
interaction with a pronounced peak in the density of state 
near the band edge |Q; 2). the Hund's rule coupling in 
the presence of the orbital degeneracy of itinerant elec- 
trons leads to ferromagnetism ||; and 3). the motion 
of electrons in the localized spin background forces the 
spins to align parallelled, i.e., the double exchange mech- 
anism j3j. Significant progresses have been made in our 
understanding of ferromagnetism in last several years. 
Recent reviews on this subject are seen in Refs. @-§ 
and references therein. On the other hand, ferromag- 
netism in Ri^jXjMnOj and related ordering states have 
stimulated extensive interests due to the phenomenon 
of the colossal magnetoresistance. The physical origins 
are still lack of full understanding. Electronic models 
for this family of materials contain all the three factors 
which favor to ferromagnetism, and are ideal candidates 
to test all physical mechanisms for ferromagnetism. Var- 
ious techniques are applied to investigate magnetism in 
these models 0. Due to the complexity of the systems 
quite few rigorous solution or rigorous results are ob- 
tained to the models. Usually rigorous results are very 
helpful to shape the physics of theoretical models, espe- 
cially in such a strong correlated electron system. 

In this Letter we shall present an exact solution of the 
ground state for a one-dimensional electronic model for 
transition-metal oxides in the case of strong coupling. 
The ground states are proven to be ferromagnetic by 
means of the Perron-Frobenius theorem. The physical 
mechanism and its relevance to three-dimensional cases 
are discussed. Due to the orbital degeneracy of itinerant 
electrons, the superexchange coupling can be ferromag- 
netic rather than antiferromagnetic in the one-band case. 

An electronic model Hamiltonian for the transition- 
metal oxides is defined on a discrete lattice A with N A 



site and is written as 



where 



H = Ht + H,. + H; + Hi, 



(1) 



Ht — J]] £y c i, lt a C 3,l' ,a\ 

z,*y ,7' ,(?,(?' 

J~ff'( C i,~f,<J Ci ,l^' C l,j',cr' C j,l',°' 
~^ C l,f,cr C i,l' ,f c I,7,er' c «,7' ,"')'■> 

Hh = — ^ JhSi ■ Sj )7 . 



H, 



H 



z,7 



and 



are creation and annihilation opera- 



z,7,<r 

tors for e g electron with spin a (= ±1) on orbital 
7(= ±1) at site i, respectively. Si is the maximal to- 
tal spin of three t 2g electrons (S = 3/2), and Si j7 = 

Zw ^CT' c l 7 ,cr c i,7,<T'/ 2 is a s P m operator for an e g elec- 
tron and a are the Pauli matrices. H t describes the pro- 
cess of electrons hopping between nearest neighbor sites. 
H u is of the the on-site Coulomb interaction. Hj is the 
Hund exchange coupling between itinerant electrons on 
different orbital at the same site and J 77 < is always posi- 
tive. Hh is the Hund coupling between localized spin and 
itinerant electrons at the same sites. In this model the 
total spin operator, 



Sfot — Si + Si 



z,7 



commutes with the Hamiltonian, [S to t,H] = 0. Hence 
the total spin and its z-component are good quantum 
numbers. The maximal total spin is S rnax — N A S + N e /2 
(N e is the number of electrons and we just consider N e < 
2iV A ). We call the state with S max a ferromagnet. 

This model has been investigated by many authors . 
To the best of our knowledge, there has been no exact 
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solution for this model until now. In order to establishes 
some rigorous results for this model we first consider a 
one-dimensional open chain and focus our attention on 
the case that (i). t[- = —tS ln >Sij±i and t > 0; (ii). 
f7 77 / = U = +00 if 7 = 7'; (iii). U 7 y = U' if 7 5^ 7'; 
(iv). J 7 y = J > for 7 ^ 7'; (v). J H > 0. It is 
worth stressing that the condition (ii) excludes the double 
occupancy of electrons on the same orbital at the same 
site, which is essential to solve the one-dimensional model 
exactly. 

To solve the model rigorously, we first investigate 
the ground state properties by means of the Perron- 
Frobenius theorem The ground state is proven to 
be non-degenerate and possesses the maximal total spin. 
In this case we can write down the ground state wave- 
function explicitly by utilizing the Bethe ansatz. The 
solution is valid for any density of electrons. Perron- 
Frobenius Theorem (for a real symmetric matrix) states: 

Let M = {rriij} be a real, symmetric, and square ma- 
trix. If the matrix M satisfies the conditions: (i). all its 
off- diagonal matrix elements are non-positive, mij < 0; 
(ii). any i and j are connected by the matrix, i.e., we can 
always find an integer n such that (M n )ij ^ 0, then the 
lowest eigenvalue is nondegenerate and all elements of 
the corresponding eigenvector can be taken to be positive. 

This theorem was successfully applied to prove the 
existence of ferromagnetism in the quantum double ex- 
change model || , the one-band Hubbard model and 
the orbitally degenerate Hubbard model pH| . In order 
to apply this theorem to the present model (Eq.(l)), we 
have to choose a suitable set of basis to write the Hamilto- 
nian in the form of a real, symmetric and square matrix 
which obeys the two conditions. Consider the system 
contains N e electrons. As the z-component of total spin, 
5 t z ot , is a good quantum number, the Hilbert space can 
be decomposed by S£ ot . Assume Ni electrons on 7 = — 1 
orbit with spin o\ < 02 < ■ ■ . < erjvi are located at sites 
X\ < X2 < ■ ■ ■ < xn x , and N2 electrons on 7 = 1 orbit 
with spin ctjvi+i < 0jVi+2 < ... < ctjvi+jv 2 are located 
at sites xat 1+ i < xn x +2 < . . . < xn c - x m — 1, . . . , iV A 
and rrii = 0, ... , 25. As Ni — N 2 is also a good quantum 
number, thus we choose specific 5 t z ot , N\ and N 2 . (We 
shall show the lowest energy state is located at Ni = N2 
if N e is even, or N\ — N2 ± 1 if N e is odd.) Denote a be 
one of the configurations of electrons and localized spins, 

\ a ) ~ C x%,— ' ' ' C x Nl ,—l,(7 t f 1 C x Nl + 1 ,l,a Nl + 1 ■ ■ ■ C X Nc ,l,CTJV e 

x (S+) mi (S+) m2 ...(S+J mN -|0) 

where the state |0) has the properties: 

!)■ c x m ,~f m .a m |0) = for any m = 1, . . . , N e ; 

2). S-\Q) = for any i = 1, . . .,N A . 

As we just construct a basis within the subspace 5 t z oi , it 

is required that 



^o t = E^+E(^- 5 )- ( 2 ) 

m— 1 i— 1 

A complete set of basis consists of all possible configu- 
rations of {{x m }, {jm}, {a m }, {nii}} with the condition 
. a in I a) represents one of the possible configurations. 
On the basis we come to show that the Hamiltonian 
satisfies the condition of non-positivity and connectivity. 

(1) . Non-positive off-diagonal elements: On the ba- 
sis we choose, the Hamiltonian can be expressed in the 
form of square matrix. All non-zero off-diagonal elements 
of the matrix (a\H\a') are: (a) —t if x n = x' n ± 1 and 
all other indices are the same; (b) —J if a n = io~' m , 
x n = x' m , 7„ = — Y m and all other indices are the same; 
(c) -Jff if m, + a n = m! i + a' m , j n = j' m , x n = x' m = i 
and all other indices are the same. All the non-zero ele- 
ments are negative if t, J, Jh > 0. Therefore the Hamil- 
tonian matrix satisfies the first condition of the Perron- 
Frobenius theorem. 

(2) . Connectivity: all basis are connected through H: 
(a), the hopping terms connect all lattice sites within 
the orbit 7; (b). the Hund coupling J connects the two 
orbits at the same site; (c). the Hund coupling Jh con- 
nects the orbits of itinerant electrons and the localized 
spins. Combination of (a), (b) and (c) shows all basis are 
connected by H. 

On the basis, the lowest energy state with 5 t z ot is ex- 
pressed in the form, 

I*) = E /({ X ™}< frm}, { ff ™}< (3) 
a 

According to the Perron- Frobenius theorem, we conclude 
that the lowest energy state with 5 t z ot is non-degenerate 
and all coefficients / can be chosen to be positive, / > 0. 

Since the lowest energy state with S% ot (Eq.(||)) is non- 
degenerate and the total spin is a good quantum number, 
therefore the state must be an eigenstate of the total spin, 

(S tot ) 2 |*> = S tot (S tot + l)|tf>. (4) 

To determine Stot, we first construct an eigenstate with 
the maximal total spin S ma x — N A S + N e /2, 

\<l>(S* ot = -N A S-N e /2)} = 

El c ii,-i,-i ■ ■ ■ c x Nl ,-i,-i c x Nl+1 ,i,-i ■ ■ ■ c ijv e ,i,-ii^0' 

{x m } 

The other 2S max eigenstates with total spin Synax and 
different z-components are expressed by utilizing spin 
SU(2) symmetry, 

|$(5 z ot )) = (S+ t ) M mS! ot = -N A S - N e /2)) 

= Ei Q ) ( 5 ) 

a 

where M = 5f ot + N A S + N e /2 and the summation runs 
over all possible configurations in the subspace of 5 t z ot . 
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The main feature of the state is that all coefficients on 
the basis in the subspace S% ot are equal. This state is not 
orthogonal to the lowest energy state |\&) if they have the 
same z-component of total spin, 

(*|<&> ^0, (6) 

since all coefficients in |W) are positive. As both states 
are eigenstates of the total spin Stot, we have 

(*|(S tot ) 2 |$) = S tot (S tot + l)(*\$) 

+ (7) 

Hence we conclude that Stot — S max , i.e., the lowest 
energy state possesses the maximal total spin. In other 
word, the state is fully ferromagnetic. 

The lowest energy state with S^ ot is fully ferromag- 
netic. According to spin SU(2) symmetry of the model, 
the ground state is (2S max + l)-fold degenerate. This 
property makes it possible to write down the ground state 
wavefunction explicitly. The coefficients in the ground 
state are independent of spin indices, i.e., 

f({x m }, {7m}, Wm}, {mi}) = g({x m }, {7 m }). (8) 

Hence the Schrodinger equation, H\*f?) = E g \^f), is re- 
duced to 

-t 9{{xx,... },{jm}) 

m,S=±l 

+ U <*ff C 1 ~ *Tn,7m)ff({ a; m},{7m}) 

n<m 

= (E g - N e J/2 - N A s)g{{x m }, { 7m }). (9) 

where E g is ground state energy of the model. In the case 
of ferromagnetism, the Schrodinger equation is reduced 
to a one-band Hubbard model with the on-site Coulomb 
interaction, f/ e // — U' — J if we use the orbit indices 
instead of usual spin indices. This equation can be solved 
exactly by means of the Bethe ansatz |l2| . The solution 
to g is expressed as 

N e 

9({xm}, {7m}) = ^[Q,P]exp[i Y k P m x Q m }> ( 10 ) 

P m=l 

where P and Q are two permutations of (1,2, ... , N e ). 
The coefficients [Q,P] have the relation, 

[Q,P]=Y%+ 1 [Q,I*]; (11) 
Y i,i+i = (sinfc„ - sinfc m )P^ +1 - iU ef f/2 

(smk n -sink m )+iU ef f/2 [ ' 

where 

P= (Pi,... ^ =n,P i+1 =m,...P N ); 
P' = {P x ,...Pi = m,P' i+l =n,...Pv). 



k n (n = 1, . . . , N A ) are determined by Eqs.(|TT]) and (|T^). 
The ground state is located in the subspace Ni — N2 if 
the number of electrons is even, or Ni = N2 ± 1 if the 
number of electrons is odd. Its lowest energy is 

N c 

E g = -2t2_^cosk n --N e J -N A J H S. (13) 

n=l 

The properties of the Lieb-Wu's solution have been dis- 
cussed extensively [[l3|. Except for the filling N e = N A , 
there is no energy gap. Hence the ground state is a metal- 
lic ferromagnet. 

We have obtained the exact solution of the ground 
state of the electronic model for transition metal oxides 
in a one-dimensional open chain for any density of elec- 
trons in the case of U = +00. The state is proven to be 
fully ferromagnetic. Here we come to discuss the physical 
mechanism of ferromagnetism and its possible relevance 
to realistic and three-dimensional systems. The model 
we discuss contains all possible factors which favor to fer- 
romagnetism: large on-site Coulomb interaction, orbital 
degeneracy of itinerant electrons and the Hund's rule cou- 
pling. In the case of J = Jh = 0, the model has highly 
spin degenerated ground state when U = +00. Any fi- 
nite U will remove the degeneracy to form a spin singlet 
state. Thus without the Hund's rule coupling the only 
strong on-site Coulomb interaction cannot drive electrons 
to form ferromagnetism. On the contrary, the non-zero 
Hund's rule coupling, i.e., J ^ and Jh 7^ 0, will also 
remove the degeneracy to form a ferromagnet. This in- 
dicates that the interplay between the on-site Coulomb 
interaction and the Hund's rule coupling plays an essen- 
tial role in the stability of ferromagnetism in the model. 

The ferromagnetism may also be survived in the case 
of finite U and large J and Jh. To see this effect, we con- 
sider a more realistic and three-dimensional model which 
is used extensively to describe the doped lanthanum man- 
ganese oxides. The final consequences are also valid to 
the one-dimensional case. The transfer integrals in the 
model for Ri-^X^MnOs are assumed to take a Slater- 
Koster form given by the hybridization between e g orbit 
and nearest oxygen p orbit |H , 

t]/ = -tr^„-tr^„-tT^, (14) 
for Tj = Yi ± x, Yi ± y and ± z, respectively, where 

'•-(V?)."- (if)— (;$)• 

7 = ±1 represent ( J ~\ = 3z 2 — r 2 orbit and ( ^ ~\ = 

x 2 — y 2 orbit of e g electron, respectively. In reality, , the 
parameters of the model for the manganites are roughly 
estimated as U Jh , J 3> t |l5| . Hence we can obtain 
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an effective Hamiltonian at N e — N A by means of combi- 
nation of the projection technique and the perturbation 
technique p6| 



H eff ~ 



U'-i 



W 



3J 
2 



S«-S J - + (S + l/2)(S + 3/2) \ d 

(5+l/2)(5 + 3/2) J ij 

v / ^-(S + l/^ N d 

^\(S + l/2)(S + 3/2)J ij 



J H S^\(S+l/2)(S + 3/2)J 



U + J H S^\ (S+ 1/2) 2 



(15) 



where 5i is a spin operator with 5 + 1/2 as a ferromag- 
netic combination of the localized spin and itinerant elec- 
tron at the same site. P^ is the projection operator for 
the orbital occupancy 



discussion of the phase diagrams will be published else- 
where. 

Summarizing, the model on an open chain is solvable in 
the case of U = +oo. The ground state is metallic ferro- 
magnetic except for the density N e /N A = 1. The charge 
and orbital degrees of freedom are determined by a Lieb- 
Wu solution for one-dimensional Hubbard model. The 
ferromagnetism can be survived in finite U and higher 
dimensional system. Due to the orbital degeneracy of 
itinerant electrons and strong on-site repulsion U, the su- 
perexchange coupling can be ferromagnetic rather than 
antiferromagnetic in the one-band case. 

This work was supported by a CRCG research grant 
at the University of Hong Kong. 



P d - =T a (l- T a ) 



with a — x,y, z, which depends on the direction of r^— Tj. 
From the effective Hamiltonian, the first term favors to 
ferromagnetic correlation and other two terms favor to 
antiferromagnetic correlation. When U = +oo or suffi- 
ciently large, the third term is suppressed. Comparing 
the first and second terms, we find that the ferromag- 
netic coupling is always stronger than antiferromagnetic 
coupling since 



v - - 

u 2 



> 



i 2 



v + + j H s 



(16) 



and Pfj always has nonnegative eigenvalues. When JjjS 
is sufficiently large and U' — J/4 is kept unchanged, the 
first term is predominant even for a finite U. However 
when U decrease to finite, an antiferromagnetic frus- 
trations is introduced. The stability of ferromagnetism 
depends on the competition between the strong Hund 
coupling and the effect of finite U . In LaMnC>3, which 
N e = N/\, a ferromagnetic and insulating phase was ob- 
served at low temperatures fl7|| . This experimental ob- 
servation is in qualitatively agreement with our analysis. 
Oppositely, if we neglect the orbital degeneracy of e g elec- 
trons and use a one-band Kondo lattice model to describe 
the sample, we always have an antiferromagnetic phase at 
low temperatures [[l8| . It implies that the orbital degree 
of freedom is very important for us to understand the 
phase diagrams of doped manganese oxides. When the 
system is doped, i.e. N e < N A , due to the strong Hund 
coupling, Jh 3> t, the motion of itinerant electrons tends 
to force the localized spins to align parallelled. Thus the 
double exchange mechanism as well as the orbital degen- 
eracy and strong U should be also responsible for the 
metallic ferromagnetism in Ri_ I X I Mn03. The detailed 
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